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Abstract 

Using first-principles calculations, we investigate the adsorption of various gas molecules (H2, 02, H2O, NH3, NO, NO2, 
and CO) on nnonolayer l\/loS2. The most stable adsorption configuration, adsorption energy, and charge transfer are 
obtained. It is shown that all the molecules are weakly adsorbed on the monolayer !\AoS2 surface and act as charge 
acceptors for the monolayer, except NH3 which is found to be a charge donor. Furthermore, we show that charge 
transfer between the adsorbed molecule and M0S2 can be significantly modulated by a perpendicular electric field. 
Our theoretical results are consistent with the recent experiments and suggest M0S2 as a potential material for gas 
sensing application. 
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Background 

Sensing gas molecules, especially toxic gas, is critical 
in environmental pollution monitoring and agricultural 
and medical applications [1]. For this reason, sensitive 
solid-state sensors with low noise and low power con- 
sumption are highly demanded. While sensors made from 
semiconducting metal oxide nanowires [2,3], carbon nan- 
otubes [4,5], etc. have been widely studied for gas detec- 
tion for some time, graphene as a novel sensing material 
has further stimulated strong interests in the research 
community since Schedin et al. [6] demonstrated that 
a micrometer-sized graphene transistor can be used to 
detect the ultimate concentration of molecules at room 
temperature, presenting a pronounced sensitivity many 
orders of magnitude higher than that of earlier sensors. 
The graphene-based sensor is actualized by monitoring 
the change in resistivity due to the adsorption or desorp- 
tion of molecules, which act as charge acceptors or donors 
[7-9]. It is shown that sensitivity of this sensor can be 
further improved through introduction of the dopant or 
defect in graphene [10-13]. Despite these achievements, 
researchers continue to seek for novel sensitive sensors 
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similar to or even more fascinating than graphene gas 
sensors. 

Recently, two-dimensional monolayer M0S2, a kind of 
transition metal dichalcogenide, has attracted increasing 
attention because of its versatile and tunable proper- 
ties for application in transistor, flexible optoelectronic 
device, photodetector, and so on [14-19]. Unlike graphene 
which lacks a band gap and needs to be engineered to 
open the gap for practical application, pristine mono- 
layer M0S2 has a direct band gap of 1.9 eV [20] and 
can be readily used to fabricate an interband tun- 
nel field-effect transistor (FET) [21-26]. In this context, 
Radisavljevic and co-workers [21] first reported a top- 
gated FET on the basis of monolayer M0S2, which pos- 
sesses a room-temperature current on/off ratio exceeding 
10^ and mobility of 200 cm^ V"^ s"^ At the same time, 
the success of graphene-FET sensors also greatly inspires 
the intensive exploration of M0S2 as a sensing material. 
Since monolayer M0S2 holds a high surface-to-volume 
ratio comparable to graphene, a M0S2 -based gas sen- 
sor is expected to have excellent sensing performance as 
well. More recently, FET sensors made from mechani- 
cally cleaved monolayer and multilayer M0S2 have been 
demonstrated, which exhibit high sensitivity for NO gas 
with a detection limit down to 0.8 ppm [27]. The supe- 
rior sensitivity for NO2 has been observed in a flexible 
FET sensor array on a polyethylene terephthalate (PET) 
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substrate based on a M0S2 channel and reduced graphene 
oxide (rGO) electrodes [28]. Compared to the rGO-FET 
sensor, this novel sensor array displays much higher sen- 
sitivity, which can even be enhanced by up to three times 
via functionalization of M0S2 with Pt nanoparticles. 

Although the M0S2-FET sensor for nitride oxide has 
been experimentally realized, the underlying mechanisms 
regarding how NO^^ molecules interact with the M0S2 
surface and affect the electronic properties are not clear. 
Moreover, the response of M0S2 upon exposure to other 
gas molecules like H2, O2, H2O, NH3, CO, etc. remains 
to be examined either. In order to fully exploit the pos- 
sibilities of a M0S2 -based gas sensor, a systematic study 
on the adsorption of gas molecules on a M0S2 surface is 
thus desired from a theoretical point of view. In this work, 
using first-principles calculations, we first determine the 
most stable configuration for gas molecules adsorbed on 
monolayer M0S2, as well as the corresponding charge 
transfer between them. Modification of the electronic 
properties of host monolayer M0S2 due to the molecule 
adsorption is then examined. Furthermore, the effect of 
an external electric field on the charge transfer is also dis- 
cussed. To the best of our knowledge, no prior theoretical 
work has been conducted on these issues. 

Methods 

First-principles calculations are performed using the 
Vienna ab initio simulation package ( VASP) [29,30] on the 
basis of density functional theory (DFT). The exchange- 
correlation interaction is treated by local spin density 
approximation (LSDA). Spin-polarized calculations are 
also carried out with generalized gradient approximation 
(GGA) in some specific cases. A cutoff energy of 400 eV 
for the plane-wave basis set and a Monkhorst-Pack mesh 
[31] of 5 X 5 X 1 for the Brillouin zone integration are 
employed. In order to eliminate the interaction between 
two adjacent monolayer M0S2, a vacuum layer larger than 
15 A is adopted in the calculations. All the structures are 
fully relaxed by using the conjugate gradient method until 
the maximum Hellmann-Feynman forces acting on each 
atom is less than 0.02 eV/A. By means of Bader analysis 
[32], charge transfer between the monolayer substrate and 
the adsorbate is obtained. The electric field in VASP is 
actualized by adding an artificial dipole sheet at the center 
of the simulation cell. 

Results and discussion 

We consider the absorption of H2, O2, H2O, NH3, NO, 
NO2, and CO on two-dimensional monolayer M0S2. A 
4x4 supercell of monolayer M0S2, with a single gas 
molecule adsorbed to it, is chosen as the computa- 
tional model. The optimized lattice constant of mono- 
layer M0S2 is 3.12 A, and consequently, the distance 
between two neighboring gas molecules is larger than 



12 A. The monolayer M0S2 consists of a monatomic Mo- 
layer between two monatomic S -layers like a sandwich 
structure, in which Mo and S atoms are alternately located 
at the corners of a hexagon. In order to determine the 
favorable adsorption configuration, four adsorption sites 
are considered, namely, H site (on top of a hexagon), Tm 
(on top of a Mo atom), Ts (on top of a S atom), and B 
site (on top of a Mo-S bond). The gas molecule is initially 
placed with its center of mass exactly located at these sites. 
For each site, configurations with different molecular ori- 
entations are then examined. Take NO as an example, 
three initial molecular orientations are involved, one with 
NO axis parallel to the monolayer and two with NO axis 
perpendicular to it, with O atom above N atom and O 
atom below N atom [see Additional file 1 for more detailed 
adsorption configurations] . The adsorption energy is cal- 
culated as E(i = -£MoS2+molecule ~ (-^MoS2 ~l~ ^molecule)? 

where £'moS2 + molecule is the total energy of M0S2 with an 
absorbed molecule and ^moS2 ^^i^^ ^molecule are the total 
energies of pristine M0S2 and isolated molecule, respec- 
tively. A negative value of Ea indicates that the adsorption 
is exothermic. 

Table 1 summarizes the calculated values of equilib- 
rium height, adsorption energy, and charge transfer for 
the adsorption of gas molecules on monolayer M0S2. 
The values for each adsorbate correspond to its favor- 
able adsorption configurations obtained at different sites. 
The equilibrium height is defined as the vertical distance 
between the center of mass of the molecule and the top 
S -layer of the M0S2 sheet. Note that the adsorption ener- 
gies are often overestimated at the LDA level, but this is 
not very essential here because we are primarily interested 
in the relative values of adsorption energies for different 
configurations and finding the most favorable one among 
them. From Table 1, we see that for both H2 and O2, the 
Tm site is found to be their most favorable site with the 
adsorption energies of —82 and —116 meV, respectively. 
The corresponding structures are shown in Figure la,b. 
Nevertheless, it seems that the two molecules adopt dis- 
tinct orientations. While H2 has an axis perpendicular to 
the monolayer, that of O2 is nearly parallel to the mono- 
layer with its center of mass on top of the Tm- H2O, NH3, 
and NO2 are preferably adsorbed at the H site, resulting 
in the adsorption energies of —234, —250, and —276 meV, 
respectively. Structures for the three systems are shown in 
Figure lc,d,f. Contrary to the configuration for H2O where 
H-O bonds adopt tilted orientation with H atoms point- 
ing at the monolayer, all the H atoms of NH3 point away 
from the monolayer. NO2 is bonded with O atoms close to 
M0S2. In our calculations, H2, O2, H2O, and NH3 fail to 
have stable configuration at the B site; this is because they 
tend to migrate to other sites during structural relaxations. 
In contrast, the configuration with the center of mass 
located at the B site is found to be the most favorable 
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Table 1 Results for gas molecules on monolayer M0S2 calculated by LDA functional 



Gas 


h 


Hsite 

Err 


AO 


h 


Tm site 


AO 


h 


Ts site 

E^ 


AO 


h 


Bsite 

E^ 


AO 


H2 


2.62 


-70 


0.004 


2.61 


-82 


0.004 


3.02 


-49 


0.008 








O2 


2.79 


-106 


0.034 


2.71 


-116 


0.041 


3.19 


-64 


0.020 








H2O 


2.59 


-234 


0.012 


2.67 


-222 


0.016 


3.13 


-110 


0.009 








NH3 


2.46 


-250 


-0.069 


2.61 


-222 


-0.051 


3.21 


-100 


-0.024 








NO 


2.68 


-195 


0.011 


2.90 


-185 


0.011 


2.88 


-152 


0.039 


2.83 


-211 


0.022 


NO2 


2.65 


-276 


0.100 








2.71 


-249 


0.119 


2.62 


-249 


0.114 


CO 


2.95 


-128 


0.020 


3.22 


-124 


0.006 


3.28 


-86 


0.016 


3.15 


-128 


0.013 



Equilibrium height between the center of mass of the molecule and the top S-layer of the M0S2 sheet {h, in A), adsorption energy {Ea, in meV), and charge transfer 
from M0S2 to the molecule (AQ, in e). Negative AQ means charge transfer from the molecule to M0S2. 



one for NO, as shown in Figure le. The corresponding 
adsorption energy is determined to be —211 meV. The CO 
molecule somewhat favors both H and B sites, giving an 
identical absorption energy of —128 meV (see Figure Ig). 
For simplicity, the configuration at the H site is chosen 
as the representative for CO. All of the following results 



for these adsorbates are obtained based on their most 
favorable configurations if not specified. 

Additionally, calculations of the gas adsorption are also 
performed using GGA functional. Different from LDA 
functional which overestimates the adsorption energy, 
GGA functional usually has a tendency to underestimate 
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Figure 1 Adsorption configurations. Top and side views of the most favorable configurations for (a) H2, (b) O2, (c) H2O, (d) NH3, (e) NO, (f) NO2, 
and (g) CO on monolayer M0S2. The blue and yellow balls represent Mo and S atoms, whereas the cyanine, red, gray, and black balls represent H, 0, 
N, and C atoms, respectively. 
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Table 2 Results for gas molecules on monolayer M0S2 calculated by PW91 and PBE functionals 

LDA GGA-PW91 GGA-PBE 
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Equilibrium height between the center of mass of the molecule and the top S-layer of the M0S2 sheet {h, in A) and adsorption energy [Ea, in meV). Positive Ea indicates 
endothermic adsorption. The LDA-calculated results are also listed for reference. 



it. Upon the application of the two kinds of function- 
als, the upper and lower bounds for adsorption energy 
and other structural properties can be obtained [8]. The 
calculated values of equilibrium height and adsorption 
energy for gas molecules on M0S2 are listed in Table 2. 
Herein, two GGA functionals, PW91 and PBE, are used 
for the purpose of comparison. Both PW91 and PBE give a 
smaller adsorption energy compared to the LDA, whereas 
they show the molecules binding at an equilibrium height 
larger than that for LDA. For most molecules (with the 
exception of NO), it seems that PW91 gives more stable 
results than PBE, with their adsorption energy difference 
approximately between —7 and —28 meV. 

Next, Bader analysis is performed to predict the charge 
transfer value. It is found that most molecules studied 
except NH3 are charge acceptors with 0.004 ~ O.le 
obtained from monolayer M0S2, whereas NH3 behaves as 
a charge donor, providing 0.069e to the monolayer. The 
charge transfer values for O2 and H2O are in good agree- 
ment with recently reported values (approximately 0.04e 
for O2 and O.Ole for H2O) by Tongay et al. [33]. Note that 
our results are somewhat similar to the previous reports 
on the adsorption of gas molecules on graphene [7] 
and carbon nanotube [34], where the gas molecules also 
behave as either charge acceptors or donors. We need to 
point out that although different methods besides Bader 
analysis may give rise to different values in determining 
the electronic charge transfer, the direction and order of 
magnitude should be the same. The mechanism of the 
M0S2-FET gas sensor for NO [27] can then be under- 
stood. Before NO adsorption, the mechanically cleaved 
M0S2 channel is an n-type semiconductor in the experi- 
ment, implying that some electrons have already existed 
in the conduction band. After NO adsorption, electron 
charge is transferred to the NO molecule, inducing a p- 
doping effect on the M0S2 channel. As a result, the chan- 
nel resistance increases and current decreases. The similar 
behavior, which has been previously reported for M0S2- 
FET devices in an O2 environment [35,36], is probably 



due to the adsorption of O2 on the M0S2 surface, which 
traps electrons and sequentially reduces the current of the 
M0S2-FET. 

To further gain insight into the molecule-monolayer 
interaction, we calculate the adsorption energy curves for 
all the studied gas molecules, wherein the height between 
the center of mass of the molecule and the top S-layer of 
the M0S2 sheet is varied between 1.5 and 5.0 A. The cor- 
responding results are given in Figure 2. It is shown that 
the curve for NO2 gives the largest adsorption energy at 
the minimum, which is three times higher than that of the 
H2 curve. At equilibrium, NH3 has a minimum height of 
2.46 A with respect to monolayer M0S2, whereas CO has 
a maximum molecule-monolayer height of 2.95 A. All the 
curves nearly reach the asymptotic value at 5.0 A. Due to 
the small adsorption energy and large separation height, 
the interaction between the gas molecules and the M0S2 
surface can thus be characterized as physisorption. 

Figure 3 presents the charge density difference images 
for these molecule-monolayer systems, calculated by the 

formula Ap = PMoS2+molecule - (PM0S2 + Pmolecule)> 




Height (Angstrom) 



Figure 2 Adsorption energy versus height. Adsorption energy 
versus height between the center of mass of the molecule and the top 
S-layer of monolayer M0S2 for all the studied molecular adsorbates. 
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Figure 3 Charge density difference plots. Charge density difference plots for (a) O2, (b) H2O, (c) NH3, (d) NO, (e) NO2, and (f) CO interacting 
witli monolayer M0S2. The red (green) distribution corresponds to charge accumulation (depletion). The isosurface is taken as 5 x 1 e/A^. 
The direction and value of charge transfer are also denoted. 
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Figure 4 Band structures. Band structures of (a) pristine, (b) 
NO-adsorbed, and (c) N02-adsorbed monolayer MoS2.The black 
(red) line corresponds to the up-spin (down-spin) bands, whereas the 
dashed green line denotes the Fermi level. 
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Figure 5 Electric field effect, (a) Representation of the applied 
perpendicular electric field, where the arrows denote its positive 
direction, (b) Variation of charge transfer as a function of electric field 
strength for NO, and NO2, adsorbed on monolayer M0S2. 
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where PM0S2 + molecule > PM0S2' and Pmolecule are the charge 
density of the molecule-adsorbed M0S2, pristine M0S2, 
and isolated molecule, respectively. The red region shows 
the charge accumulation, while the green region repre- 
sents the charge depletion. It is shown that the M0S2 
sheet is considerably polarized upon the adsorption of 
gas molecules, and electrostatic interaction plays a role 
in the attractive interaction. The polarization in the H2O, 
NH3, NO, and NO2 cases are stronger than that in the 
O2 and CO cases, giving rise to a larger interaction 
energy. It explains why the former gives larger adsorp- 
tion energies (-234, -250, -211, and -276 meV for H2O, 
NH3, NO, and NO2, respectively) than the latter (-116 
and —128 meV for O2 and CO, respectively) mentioned 
above. 

We examine the electronic properties of monolayer 
M0S2 adsorbed with gas molecules. The band structure 
before adsorption is presented in Figure 4a. It is found that 
the pristine monolayer M0S2 is a semiconductor with a 
direct band gap of 1.86 eV at K point, which is in good 
agreement with reported works [37-39]. The band struc- 
tures for both valence bands and conduction bands of 
monolayer M0S2 are not significantly altered when H2O, 
NH3, and CO are adsorbed, and the gap values remain 
around 1.86 eV (not shown here). The situation is similar 
in the cases of O2, NO, and NO2 except the flat impurity 
states in the gap of the host monolayer induced by these 
adsorb ates. While O2 introduces two close-lying down- 
spin states 0.519 and 0.526 eV above the Fermi level (Ep) in 
the band gap, NO2 introduces an unoccupied down-spin 
state 0.31 eV above Ep, as given in Figure 4c. Three impu- 
rity states emerge inside the band gap upon the adsorption 
of NO, namely, one occupied up-spin state 0.12 eV below 
Ep, one unoccupied up-spin state 0.11 eV above Ep, and 
one unoccupied down-spin state close to the conduction 
band edge with an energy separation of 0.064 eV between 
them (see Figure 4b). The adsorption of O2, NO, and 
NO2 on the M0S2 surface, on the other hand, creates 
magnetic moments of 2.0, 1.0, and 1.0 /xb per supercell, 
respectively. 

As the charge transfer between the adsorbed molecule 
and monolayer M0S2 plays a crucial role in determining 
the performance of the M0S2 sensor, it may be sensi- 
tive to the applied electric field, similar to the case of 
graphene [40] . For brevity, NO and NO2 adsorbed mono- 
layers are chosen as the representative systems. Figure 5a 
gives the schematic illustration of the electric field applied 
in our study, which is perpendicular to the plane of mono- 
layer M0S2 with its positive direction aligned upward. The 
variation of charge transfer with respect to the electric 
field is shown in Figure 5b. It is found that the charge 
transfer from the monolayer to the adsorbed molecule 
increases with the increment of field strength along a pos- 
itive direction, whereas it tends to decrease once reverse 



electric field is applied. This negative electric field will 
drive the electrons from the molecule to the monolayer 
when its field strength is beyond a critical value. While 
NO and NO2 attain 0.022e and O.le in the absence of elec- 
tric field (£ = 0 V/A), respectively, they turn out to accept 
0.225e and 039e from monolayer M0S2 at £ = 1 V/A 
and conversely donate 021e and 0.028e at £ = —1 V/A. 
The dependence of charge transfer on field direction is 
probably attributed to the dipole moment of the molecule- 
monolayer system [41]. Band structure calculations for 
the two systems, on the other hand, show that the impu- 
rity states in the band gap shift towards the valence or 
conduction bands of monolayer M0S2, depending on the 
direction of the applied perpendicular electric field. 

Conclusions 

In this work, we present a first-principles study on the 
structural and electronic properties of monolayer M0S2 
upon adsorption of gas molecules. Various adsorption 
sites and molecule orientations are involved to determine 
the most stable configurations. We find that all molecules 
are physisorbed on monolayer M0S2 with small charge 
transfer, acting as either charge acceptors or donors. Band 
structure calculations indicate that the valence and con- 
duction bands of monolayer M0S2 is not significantly 
altered upon the molecule adsorption, though certain 
molecules such as O2, NO, and NO2 introduce adsor- 
bate states in the band gap of the host monolayer. In 
addition, we demonstrate that the application of a perpen- 
dicular electric field can consistently modify the charge 
transfer between the adsorbed molecule and the M0S2 
substrate. Our theoretical findings show that M0S2 holds 
great promise for fabricating gas sensors. 
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